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Abstract 



We study a class of oscillating bounce solutions to the Euclidean field equa- 
tions for gravity coupled to a scalar field theory with two, possibly degenerate, 
vacua. In these solutions the scalar field crosses the top of the potential bar- 
rier k > 1 times. Using analytic and numerical methods, we examine how 
the maximum allowed value of k depends on the parameters of the theory. 
For a wide class of potentials fc max is determined by the value of the second 
derivative of the scalar field potential at the top of the barrier. However, in 
other cases, such as potentials with relatively flat barriers, the determining 
parameter appears instead to be the value of this second derivative averaged 
over the width of the barrier. As a byproduct, we gain additional insight 
into the conditions under which a Coleman-De Luccia bounce exists. We dis- 
cuss the physical interpretation of these solutions and their implications for 
vacuum tunneling transitions in de Sitter spacetime. 
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I. INTRODUCTION 



There has recently been renewed interest in the problem of vacuum tunneling in de Sitter 
spacetime. Although this interest has been sparked in large part by developments in string 
theory, the tunneling problem itself can be addressed within the context of quantum field 
theory. The essence of the problem is captured by considering a theory with a single scalar 
field described by the Lagrangian 

c= l -{d,ct>f-v{ ( t>) (i.i) 

where the scalar field potential has two unequal minima, as illustrated in Fig. I. The 
lower minimum corresponds to the stable "true vacuum" state, while the higher minimum 
corresponds to a metastable "false vacuum" . 

At zero temperature, and in the absence of gravitational effects, the false vacuum decays 
via a quantum mechanical tunneling process that leads to the nucleation of bubbles of true 
vacuum. The semiclassical calculation of the bubble nucleation rate per unit volume, T, is 
well understood [1,2]. It can be written in the form 

T = Ae~ B (1.2) 

where B is obtained from the action of the "bounce" solution to the Euclideanized field 
equations. This bounce solution has a region of approximate true vacuum (essentially, a 
four- dimensional bubble) separated by a wall region from a false vacuum exterior. 




FIG. 1. The potential for a typical theory with a false vacuum. 

At finite temperature, bubble nucleation proceeds via both quantum tunneling and ther- 
mal fluctuations. When the temperature T is high enough that the latter process dominates, 
T can still be written in the form of Eq. (1.2), but with B = E/T, where E is the energy of 
a critical bubble. 

In this paper we will be concerned with vacuum tunneling in situations where gravita- 
tional effects are important. We will assume that V(<j>) > everywhere, so that both the 
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false and true vacua correspond to de Sitter spacetimes. This problem was first discussed 
by Coleman and De Luccia [3] , who argued that the proper generalization of the flat space- 
time calculation of B could be obtained by looking for a bounce solution of the Euclidean 
version of the coupled matter plus gravity field equations. When the relevant mass scales 
are well below the Planck mass and the bubble sizes are small compared to the curvature 
of the de Sitter space, the Coleman-De Luccia prescription leads to a description of bubble 
nucleation that is very similar to that of the zero-temperature flat spacetime case, with small 
gravitational corrections of the expected order of magnitude. 

However, when the mass scales are higher or the bubble sizes larger, not only are the 
quantitative deviations from the flat space case more significant, but there are also qualita- 
tive differences that raise issues of interpretation. These suggest that the de Sitter vacuum 
transition process has aspects of both quantum tunneling and of thermal fluctuation. Fur- 
ther, the nucleation of true vacuum bubbles in false vacuum regions seems to imply the 
possibility of nucleating false vacuum bubbles in true vacuum regions, with the relative 
rates of the two processes having a natural thermal interpretation [4] . 

One of the most striking consequences of the inclusion of gravity is the existence of a 
homogeneous Euclidean configuration, the Hawking-Moss solution [5], that is quite different 
in form and physical interpretation from the Coleman-De Luccia bounce, and that has no 
counterpart in the flat spacetime problem. 

In this article we will focus on yet another type of Euclidean solution, which might be 
termed an oscillating bounce. In contrast with the Coleman-De Luccia bounce, where the 
scalar field varies monotonically from 0f v toward tv , and the Hawking-Moss solution, where 
it is everywhere equal to its value at the top of the barrier, 0top, the field in these solutions 
oscillates back and forth between the two sides of the potential barrier. Like the Hawking- 
Moss solution, these have no finite action counterparts in flat spacetime. Solutions of this 
type have also been discussed by Banks [6]. Here we consider these in more detail, and 
examine the conditions under which they can exist. We find (in disagreement with [6]) that, 
for fixed values of the parameters of the theory, there are only a finite number of oscillating 
bounce solutions. 

In Sec. II, we review the formalism for calculating V in flat space-time, both at T = and 
at finite temperature. We emphasize in particular the aspects that elucidate the physical 
meaning of the bounce solution, and discuss the path integral derivation of the formula for 
the prefactor A in Eq. (1.2). In Sec. Ill we review how the formalism is generalized to include 
gravitational effects, and outline the major features of the Coleman-De Luccia analysis. With 
this preparatory material behind us, we begin our discussion of new solutions in Sec. IV. We 
develop a framework for the analysis, and then first apply it to "small amplitude" oscillating 
bounces, for which a linearized approximation can be used. A critical role is played here by 
a parameter (3 that measures the second derivative of the potential at the top of the barrier 
relative to the curvature scale of the de Sitter spacetime. For a broad class of models, the 
number of such small amplitude solutions increases with (3, with new bounces appearing as 
(3 increases through certain critical values, and with no such solution at all if (3 is too small. 
The bound obtained here is physically plausible and consistent with previous discussions; 
however, as we will see, it is not universally applicable. In Sec. V we draw some intuition 
from the analysis of the previous section and use it to discuss bounce solutions, lying outside 
the small amplitude regime, for which the effects of the nonlinear terms are dominant. We 
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then test this intuition by numerical solution of the bounce equations. In Sec. VI we turn 
to the case potentials that are unusually flat at the top, for which the relation between /3 
and the number of bounce solutions that was found in Sec. IV does not apply. In the course 
of studying these we will gain some insight into the conditions under which a Coleman-De 
Luccia bounce can exist. We will show by explicit example that bounce solutions can exist 
even if the potential is absolutely flat at the top of the barrier, despite some suggestions to 
the contrary in the literature. Finally, in Sec. VII we discuss the physical interpretation of 
the oscillation bounce solutions and include some concluding remarks. 



II. REVIEW OF VACUUM DECAY IN FLAT SPACETIME 

A. Flat space tunneling at zero temperature 

We first recall some results concerning quantum tunneling in a system with more than one 
degree of freedom. Consider a system with coordinates qi (j = 1, . . . , N) whose dynamics 
is governed by the Lagrangian 




(2.1) 

and let the point q^ v be a local (but not global) minimum of the potential energy. Given 
a system whose wave function is initially localized about qt v , we want to know the rate at 
which the system tunnels through the surrounding potential barrier. 

For each path q(s) through the potential barrier that begins at q(si) = qt v and emerges 
from the barrier at a point q(s2) = q2, one can calculate a one-dimensional WKB tunneling 
integral 



fS-2 

B 

Isi 



2 r dsyj2m[U(q(s)) - U iv ] , (2.2) 

J Si 



where £/f v = U{q^) = U{qo). The leading WKB approximation to the tunneling rate is 
proportional to e~ B , evaluated along the path that minimizes the tunneling integral [7]. 
The end point q 2 of this path is the most probable place for the system to emerge from the 
barrier, and thus gives the initial condition for the classical evolution of the system after 
tunneling. 

By manipulations analogous to ones familiar from classical mechanics, the problem of 
minimizing the integral in Eq. (2.2) can be recast as the problem of finding a stationary 
point of the Euclidean action 



T2 

dr 

Tl 



(2.3) 



The solution %(r) of the Euclidean Euler-Lagrange equations gives the same path as before, 
although with a different parameterization. Because q^ v is a minimum of U, the initial 
Euclidean time must be taken as t\ = — oo; the final time r 2 is arbitrary. 
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Since dq^/dr vanishes at T2, this solution can be continued, in a "r-reversed" form, to 
give a "bounce" solution that runs from q {v to q 2 and then back to q^ v . The exponent in the 
tunneling factor is then 



B 



dr 



m ( dql 



2 \dr 

S E (qb) - S E (qtv) 



+ U(q h ) - U{q h ) 



(2.4) 



where the factor of 2 in Eq. (2.2) has been absorbed by the doubling of the action that 
results from considering the full bounce. 

Adapting [1] this formalism to the scalar field theory of Eq. (1.1), one is led to consider 
the Euclidean action 



S E = J dx 4 d 3 x i ^ 



J 



d x 



-{d a <pf + v{<p) 



and to seek a bounce solution of the Euclidean field equation 

O a da<t> = TT ■ 



(2.5) 



(2.6) 



To match the initial state before tunneling, this solution must tend to 0f v as x 4 — > ±00, 
while the finiteness of the tunneling exponent 

B = S E ((f>bouncc) - S E (<j>{ v ) (2.7) 

requires that it also tend to 0f v at spatial infinity. The bounce contains an interior region, 
in which the field is on the true vacuum side of the barrier, that is separated from the false 
vacuum exterior by a wall of finite thickness. A spatial slice through the center of this 
solution gives the three-dimensional configuration that is both the optimal end point of the 
quantum tunneling and the initial condition for the subsequent classical evolution. This 
configuration contains a true vacuum bubble embedded in the false vacuum, with the total 
potential energy (i.e., the sum of the gradient energy and the scalar field potential) being 
equal to the initial false vacuum energy. 

This WKB calculation only gives the exponential factor in the tunneling rate. The pre- 
exponential factor is most easily obtained by path integral methods [2]. The basic strategy 
is to use a Euclidean path integral to calculate the energy of the false vacuum state. Since 
this is an unstable state, its "energy" is complex, with its imaginary part giving the decay 
rate. 

Specifically, consider the path integral 



j(t) = / m 



-S E (4>) 



(2.8) 
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where the integration is restricted to paths obeying 0(x, — T/2) = 0(x, T/2) = 0f v . As 
T — > oo the path integral is dominated by the lowest energy state with nontrivial overlap 
with these boundary conditions; i.e., by the false vacuum. Hence 

hi/(T) . . 

£ fv = -lim— ^ (2.9) 



and the bubble nucleation rate per unit volume is 

"Im ln/(T) 



T = 2 lim 



(2.10) 



where Q is the volume of space, assumed to be taken to infinity at the end of the calculation. 

The path integral can be evaluated by summing the contributions from its stationary 
points. The first of these is simply a constant homogeneous false vacuum configuration, 
<f)(x) = </>&. To leading order, the contribution of this to the path integral is 

J (T) = e~ s ^) [det S'U^T 1/2 , (2.11) 

with S denoting the second variation of the action about the (constant) classical 
solution. 

Next is the bounce solution, (f>t,(x). In calculating the determinant factor here, one finds 
that S'^(4>h) has four 1 zero modes and one negative mode. The former must be replaced 
by collective coordinates specifying the location in space and Euclidean time of the bounce, 
while the latter gives a factor of i when the square root of the determinant is taken. The 
contribution to the path integral is 

h{T) = % -QTJ |det'^(0 b )|- 1/2 e~ s '^ = % -QTK e~ B J (T) , (2.12) 

where J is the Jacobean factor from replacing the zero modes by collective coordinates, the 
factor of VtT is from the integration over the collective coordinates, the 1/2 comes from 
a careful treatment of the negative mode integration, and the prime on the determinant 
indicates that it does not include the zero modes. 

Finally, there are the contributions from the approximate stationary points corresponding 
to many widely separated bounces. 2 These are of the form 



UT) = - 



1 \i 



-VtTK 
2 



e- nB I (T), (2.13) 



where the factorial enters because interchanging bounces does not give a new contribution. 
Summing over n gives 



1 In theories with internal symmetries there may be additional zero modes; for an example of this, 
see Ref. [8]. 

2 The dilute gas approximation of considering only widely separated bounces is valid for B S> 1. 
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J(T) = £ 4(T) = 7 (T) exp 



n=0 



-ttTKe- B 



(2-14) 



The factor of io(^~) onr y contributes to the real part of the energy, so Eq. (2.10) gives 

r = fsr e - B . (2.15) 

The bounce equation is usually solved by assuming 0(4) symmetry, so that is a function 
only of s = ^x 2 + x\. Equation (2.6) then reduces to 



d 2 (j) 3 d(j) _ dV 
ds 2 s ds dd> 



(2.16) 



with the boundary conditions 

d(f> 
ds 



= , 0(oo) = fv , 



(2.17) 



s=0 



The analysis of this equation is aided by considering an analogous problem in which s 
is time and <fi represents the position of a unit mass particle with potential energy U = —V 
that is subject to a frictional force proportional to 3/s. Finding a solution is tantamount 
to finding an initial "position" 0(0) such that as s — > oo the particle comes to rest at 0f v , 
where U = —V has a local maximum. The existence of such a position can be established by 
Coleman's "overshoot/undershoot" argument [1], which relies on the fact that the friction 
decreases monotonically with time. 

Several points should be noted: 

1) The fact that the fluctuations about the bounce included only a single negative mode 
played an essential role in the path integral derivation. The false vacuum energy would have 
been real if there had been an even number of negative modes, while for an odd number of 
the form 4k + 3, the imaginary part of the energy would have had the wrong sign. 

2) Although the 0(4) symmetry of the Euclidean Lagrangian is technically quite useful, 
it does hide a very real asymmetry in the physical interpretation of the x a . While xi, x 2 , 
and X3 are ordinary spatial variables, £4 is simply a convenient parameterization of the 
configurations 0(x; x 4 ) that define the optimal tunneling path. 

3) Although the one-bounce solution has a natural interpretation in terms of the optimal 
tunneling path, there seems to be no analogous simple interpretation for the multi-bounce 
solutions. 



B. Flat space tunneling at finite temperature 

This formalism can be readily extended [9,10] to the case of finite temperature T. Instead 
of obtaining the bubble nucleation rate in terms of the imaginary part of the energy of the 
false vacuum, one instead calculates it in terms of the imaginary part of the free energy 
of the false vacuum. This can be obtained by recalling that the partition function is given 
by a path integral over configurations that are periodic in Euclidean time with period 1/T. 
Hence, the bounce solutions satisfy the Euclidean field equations on R 3 x S* 1 rather than on 
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R 4 . The only boundary conditions are at spatial infinity, where is required to take on its 
false vacuum value. 

Two quite different types of periodic bounce solutions immediately come to mind. At 
low temperatures (i.e., when 1/T is much greater than the characteristic bounce radius), a 
periodic bounce can be obtained by a small deformation of the zero temperature bounce. 
The action of this solution, which clearly corresponds to quantum tunneling, differs only 
slightly from that when T = 0. 

In the second type of bounce, corresponding to bubble nucleation by thermal fluctua- 
tions, the fields are independent of Euclidean time and so trivially satisfy the periodicity 
conditions. Since the integral over the Euclidean time is trivial, the action can be written 
as S = E/T, where E is the free energy of the static solution of the three-dimensional field 
equations. 3 These thermal bounces are subdominant at low temperature, but dominate at 
high temperature. 

It is instructive to contrast the spatial slices of the vacuum tunneling bounce with those 
of the thermal bounce. In the former case, the hypersurface at x 4 = — oo gives the initial 
false vacuum configuration before tunneling, and a slice through the center of the bounce 
gives the optimal exit point from the potential energy barrier, which is a configuration 
containing a supercritical bubble that will expand from rest. In the latter case, the spatial 
hypersurfaces give a static configuration containing a critical bubble, whose radius is at the 
top of the potential energy barrier, balanced between expansion and contraction. The initial 
state is only evident from the asymptotic behavior at spatial infinity. 



III. VACUUM TUNNELING IN CURVED SPACETIME 



Coleman and De Luccia [3] extended the bounce formalism to include the gravitational 
effects on vacuum decay. They argued that one should an Einstein-Hilbert term to the 
Euclidean action and then seek bounce solutions of the resulting Euclidean field equations. 
The tunneling exponent would then again be given by Eq. (2.7), but with the actions now 
including the additional gravitational terms. Their discussion did not, however, include the 
calculation of the prefactor A, an issue that remains poorly understood. 

If one assumes 0(4) spherical symmetry, the metric can be written in the form 



ds 2 = de + P (0 2 dn 2 3 , 



(3.1) 



where dQ 2 is the metric on the unit three-sphere, and all scalar fields depend only on £. The 
action reduces to 



Se = 2?r 2 / d£ L 3 Q0 2 + V) + ^ (p 2 p + p (? - p) 



(3.2) 



3 To be more precise, E is actually a temperature-dependent quantity E{T). In the path integral 
formalism this temperature-dependence comes about from a reordering of counterterms in the 
Lagrangian, corresponding to the standard field theory calculations [11,12] of the high temperature 
effective potential. 
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with overdots denoting differentiation with respect to £. 
The Euler-Lagrange equations are summarized by 



p aq> 



and 



Because of the second term in Eq. (3.3), must vanish at the zeros of p. If V((f>) > 0, 
which we will assume, there are always two such zeroes, so that the Euclidean manifold 
is topologically a four-sphere [13]; the only exception occurs if V{4>) vanishes at one of its 
minima, in which case the homogeneous solution with at this minimum has the flat metric 
of R 4 . One of the zeroes of p can be taken to be £ = 0. If the second zero is denoted £ max , 
then the boundary conditions on are that 

0(O) = 0(£ max ) = O. (3.5) 

The symmetry of these boundary conditions is in contrast with the flat space boundary 
conditions of Eq. (2.17). 

One consequence of the four-sphere topology is that the Euclidean actions of the bounce 
and the false vacuum are separately finite. In flat space, Se^v) is formally infinite, and 
Eq. (2.7) only makes sense in terms of a point-by-point subtraction of V((j){ v ) inside the 
spatial integral. In contrast, the Euclidean false vacuum solution of Eqs. (3.3) and (3.4) is 
a four-sphere of radius Hf 1 , with 

p(0 = Hf 1 sin(# f (3.6) 

and 



H f = 

Its action is 



(3.7) 



with the factor of —3/8 being the sum of +3/8 from the matter term in Eq. (3.2) and —3/4 
from the second, gravity, term. 

One would expect gravitational effects on tunneling to be small if the characteristic 
scales of V are much less than the Planck mass, and the flat space bounce radius f <C Hf 1 . 
In this regime, the Coleman-De Luccia bounce solution has a central true vacuum region, 
< £ ;$ f, where the scalar field profile closely approximates that of the flat space bounce. 
Outside this region <p rapidly approaches its false vacuum value although, because £ max is 
finite, 0b (0 — 0fv never quite vanishes. The tunneling exponent B differs from the flat space 
value by a fractional amount of order f 2 Hf. 
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FIG. 2. A Coleman-De Luccia bounce solution for the case where the flat space bounce radius is 
much less than Hf 1 . The picture should be visualized as a four-sphere viewed head-on, with £ = 
being the point at the center of the cross-hatched region where the field is on the true vacuum side of 
the potential barrier; £ = £ max is the antipodal point on the opposite side of the sphere. The dashed 
line passing through £ = (and also through £ = £ max ) denotes a three-sphere corresponding to 
the spatial hypersurface on which the bubble materializes. The three-sphere denoted by the lower 
dashed line is roughly analogous to the false vacuum initial-state hypersurface in the flat space 
problem. 

In this regime, the physical interpretation of the bounce solution carries over with only 
slight modifications from the flat space case. As before, a three-dimensional slice through 
the center of the bounce (see Fig. 2) can be viewed as giving initial data for the classical 
evolution; the fact that this slice has finite volume reflects the finite volume of (most) 
spacelike slices through de Sitter spacetime. In flat space, the hypersurface at £4 = —00, 
with <f>(x) identically equal to its false vacuum value, represented the initial state before 
tunneling. Although this has no precise counterpart here, one can view a slice such as that 
shown by the lower dashed line in Fig. 2 as being analogous to a slice at large negative x 4 
in the flat case. There are also multibounce configurations that are approximate stationary 
points of the action, just as in the flat space case. Summing the contributions of these to the 
Euclidean path integral, one would obtain an exponential of the one-bounce contribution 
(modulo the difficulties of properly defining the fluctuation determinant), just as in the 
previous case. 

There is, however, one difference, related to the point £ = £ max , that should be pointed 
out here. One might have expected this to be analogous to s = 00 in the flat space bounce. 
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In the latter, s = oo defines a three-sphere and corresponds both to spatial infinity on 
the hypersurface where the bubble materializes and to the initial state hypersurface. By 
contrast, £ = £ max is a single point that lies on the spacelike hypersurface that gives initial 
data for the Lorentzian problem. However, as can be seen from Fig. 2, it does not lie on the 
initial state hypersurface. 

The deviations from the flat space case become quite significant when the flat space 
bounce radius is comparable to or greater than Hj 1 . In some such cases, there is no Coleman- 
De Luccia bounce at all. When there is a bounce, its true vacuum region (whose size may 
be very different than in the flat space case) occupies a significant fraction of the Euclidean 
four-sphere, and one can no longer identify even an approximate initial-state hypersurface. 
Indeed, the asymmetry between the true vacuum interior and the false vacuum exterior is lost 
to a large degree. This suggests that the two regions can be viewed as being on equal footings, 
and that the bounce can describe either the production of a true vacuum bubble in a false 
vacuum background or that of a false vacuum bubble in a true vacuum background, with the 
initial state determined not by the properties of the bounce solution, but by whether Se^^v) 
or Se^u) is subtracted when calculating B. (A detailed examination of this possibility 
shows that the rates for nucleation of false vacuum bubbles in true vacuum and true vacuum 
bubbles in false vacuum are related by a factor that has a natural thermal interpretation, 
and that takes a simple Boltzmann form in some limits [4].) 

Finally, note that any multibounce solutions are clearly limited and quite constrained 
in this large-bounce regime, so there is no exponentiation of the bounce factor; this fact by 
itself gives a hint of the difficulties of generalizing the path integral calculation of Y. 

Whether or not there is a Coleman-De Luccia bounce solution, there is always a second 
type of solution, first pointed out by Hawking and Moss [5], that has no analogue in flat 
spacetime. Like the pure false vacuum solution, this is homogeneous, but with the scalar 
field everywhere equal to 0top, its value at the top of the barrier. Calculations exactly 
analogous to those leading to Eqs. (3.6)-(3.8) give 



p(0 = H~l sin(tf top O (3.9) 



and 



3 

Bum = S E ((j)top) - S E (<j>f v ) = g 



V(0top) Vifc) 

with if top defined by a formula analogous to Eq. (3.7). 



(3.10) 



IV. NEW SOLUTIONS 

Our focus in this paper is on a third class of Euclidean solutions that are neither simple 
bounces nor Hawking-Moss. We continue to assume 0(4) symmetry, so Eqs. (3.1)-(3.5) still 
apply. However, we now allow <fi to cross the barrier an arbitrary number of times between 
£ = and £ max ; we will denote this number by k, so that k = 1 corresponds to the Coleman- 
De Luccia bounce. The possibility of multiple barrier crossings is in contrast with the flat 
space case, where only k — 1 solutions can exist. 

To simplify matters, we will restrict ourselves to the case of a potential with only two 
local minima. Without any loss of generality we can assume that </>f v < (f) top < <p tv . 
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Our problem can be viewed as that of finding values of 0o for which the initial conditions 
0(0) = 00) 0(0) = p(0) = imply that 0(£ ma x) = 0. By scanning the range 0f v < < tv 
for such critical values, all bounce solutions can be found. In fact, because the problem 
is invariant under the transformation £ — > £ max — £, every solution with odd k will be 
encountered twice in such a scan, once as a "false-to-true" bounce and once as a "true-to- 
false" one. There may also be false-to-false and true-to-true solutions for which the initial 
and final values of are on the same side of the barrier. These appear either once or twice 
during the scan, depending on whether or not the two endpoints are identical; we have found 
examples of both types of behavior. 

It should be stressed that our use here of phrases such as "false-to-true" refers only to 
the variation of as £ ranges from to £ max - As was pointed out in the previous section, 
the value of the field at £ max [in contrast to 0(oo) in the flat space bounce] is not necessarily 
related to the initial state. 

In the flat space problem, the possible values of O can be divided into undershoot 
and overshoot regions, with the critical bounce value lying at the boundary between the 
regions. In the curved space case, where trajectories can cross top more than once, this 
overshoot /undershoot classification becomes somewhat ambiguous. Instead, we find it more 
useful to define a function h(4>o) whose zeroes correspond to the critical trajectories. 

An obvious candidate for this function would be the value of 0(£ ma x) ° n the trajectory 
specified by O . However, this diverges whenever it is nonzero, and so does not give a 
satisfactory /i(0o). The quantity p 3 is better behaved, but it too can diverge if the trajectory 
goes into a region where dV/dcj) is unbounded. However, we can avoid any such divergence 
by modifying the growth of V at large 0. This will have no effect on the critical trajectories, 
since these remain within the finite interval 0f v < < tv . 

Let us therefore define h(<f> ) to be the value of p 3 at £ max , evaluated on the trajectories 
of the modified potential. The zeros of h are the values of 0o that give bounce solutions 
satisfying the boundary conditions. Although h depends in detail upon precisely how the 
potential was modified, the locations and nature of its zeros do not, and it is really only these 
with which we are concerned. There will always be zeros at 0f v , tv , and top , corresponding 
to the pure false vacuum, pure true vacuum, and Hawking-Moss solutions, respectively. 
Depending on the potential these may be the only zeros of h, or there may be others. 

The plot of h(4>o) evolves as the parameters of the theory are varied. During the course of 
this evolution the locations of the zeros will vary continuously. The function may also gain 
or lose zeros, but always in pairs. This can happen either at top or at a point where h((f) ) 
momentarily develops a double zero. Following the motion, appearance, and disappearance 
of these zeros provides a framework for the study of the new bounce solutions. 

We will start our investigation by considering "small amplitude" solutions in which is 
confined to the region near 0t O p, where a linear approximation to Eq. (3.3) can be applied. 
We will then follow the evolution of these solutions as a change of parameters takes them 
outside the linear regime. For some potentials, all bounce solutions arise in this manner. 
Other potentials admit a second type of solution, that are never confined to the linear region. 
These can coexist with the former type of solution, or they may be the only type of bounce 
present. 

We begin by introducing an approximation that considerably simplifies the analysis of 
the bounce equations. 
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A. The fixed background approximation 



For the most part, we will focus on the case where the fractional variation of V(<f>) is 
small as varies over the range from fv to <f> tv . This allows us to write 



V(<j>) = V + V{ct>) 



(4.1) 



with |V(0)| <C Vq for all relevant values of 0. To leading approximation, the curvature of 
the Euclidean space is then independent of the value of 0, and we can write 



with 



2 _ 8tt\/ 



H Z EE 



3M 2 Pl 



(4.2) 



(4.3) 



If we define y = the equation for becomes, to leading approximation, 

d 2 <p 



, + 3 cot y 
dy 2 



dy 



1 dV 
lP~d4 



(4.4) 



with the boundary conditions that d(f>/dy vanish at both y = and y — ir. 

Once a solution to Eq. (4.4) has been found, the first correction to the metric can be 
obtained by substituting this solution into a linearized version of Eq. (3.4) and solving for 
5p = p — po- Because po is a solution of the zeroth order problem, the terms linear in Sp 
do not contribute to the action. The terms quadratic in 5p are suppressed relative to the 
matter terms in the action by a factor of order V /V , allowing us to write 



Sr — 



3 Mi 



8 Vn 



+ 2tH / d^ 



1 



+ V) + 



(4.5) 



When the false vacuum action is subtracted from that of the bounce, the first terms on the 
right hand sides cancel, leading to 



B 



9 M 4 



pi 



32 Vn 



dy sin 3 y 



dy 



(4.6) 



In particular, for the Hawking-Moss solution 



Bum — — 



pi 



3 M 4 



[v(0to P ) - v(M ■ 



(4.7) 
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B. Linearized equations 



We begin our analysis by seeking small amplitude solutions where 0o is near <p top and (j> 
remains everywhere close to the top of the barrier. These are essentially small perturbations 
about the Hawking-Moss solution. 

We expand the scalar field potential as 



V(<f>) = v{<p top ) + H 2 



-f (0 " 0to P ) 2 + \{<t> ~ 0to P ) 3 + \{<t> ~ 0to P ) 4 + 



(4.8) 



Because top is a maximum of the potential, f3 is necessarily positive, and 



(4.9) 

The sign of b can be reversed by a simple redefinition of 0, and hence is not physically 
significant. Finally, A can take either sign, although it must be positive if there are no 
higher order terms in the potential. 

If the field remains sufficiently close to top , the cubic and higher terms in the potential 
can be dropped in a first approximation, so that Eq. (4.4) takes the form 

° = + 3 COt V Ty + M ~ 0top) • (410) 



If we write <f>(y) — top = f(y) sin 3 ^ 2 y, this becomes 



Except near the outer edges of the interval < y < ir, the last term in the brackets can be 
ignored, and we find that 

4>{V) ~ 0to P « A ' L , (4.12) 

sin ' y 

with A and S constants. 

More precisely, we can note that an exact solution of Eq. (4.10) is given by a Gegenbauer, 
or ultraspherical, function 

<P(y)-<Pto P = AC 3 J 2 (cosy) (4.13) 

with a the positive root of 

a(a + 3) = (3. (4.14) 

The fact that C^/ 2 (l) is finite for all real positive a guarantees the vanishing of d<p/dy at 
y — 0; it is this condition that eliminates the Gegenbauer function of the second kind, 
Dl/ 2 (cosy). The second boundary condition, that dcp/dy also vanish at y = it, is only 
satisfied if a is an integer. Hence, the linearized approximation, Eq. (4.10), has an acceptable 
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solution only if V"((f) top ) = N(N + 3)H 2 for N = 1, 2, . . .. Because C 3 J 2 (cosy) has iV zeros in 
the range < y < it, these solutions are bounces with with k = N crossings of the barrier. 

We recall, for later use, that the Gegenbauer functions for integer N are polynomials. If 
defined with the standard normalization 

*D = ™, (4.15) 
they obey the orthogonality relation 

f 1 A n 2x^3/2,^3/2,, , 2{N+l){N + 2) _ 

y_ i dx(l-x )C M (x)Cy (x) = 5m,tv ^ N - ^ = $m,n k N . (4.16) 



C. Incorporating nonlinearities 

The linearized equations only have solutions obeying the boundary conditions when /3 is 
equal to one of a discrete set of critical values. This condition is relaxed once the nonlinear 
terms are included. This is quite similar to the way in which adding anharmonic terms 
to a simple harmonic oscillator allows small amplitude oscillations with a continuous range 
of periods. Just as the amplitude is related to the period in the oscillator problem, the 
amplitude of the bounce solutions is determined by the distance A = (3 — N(N + 3) from 
the critical value. 

Because the Gegenbauer polynomials form a complete set, an arbitrary function cf)(y) 
obeying the boundary conditions of Eq. (3.5) can be expanded as 

1 oo 

m = 0to P + ^=J2 A M C% 2 {y) , (4.17) 

A M=0 



where the coefficients Am are dimensionless. Substituting this into Eq. (4.4) and retaining 
the contributions from the cubic and quartic terms in the expansion of the potential yields 

oo ! t 

= E c3 M(y) { IP - M(M + 3)}A M - -= E A I A.j Pl j, M - sgn (A) £ AjAjA K q IJK;M 

M=0 I A / 1 A I I, J IJ,K 



(4.18) 

where the terms involving 

Pij-M = kll 12 £ dy cT{y)CT (y)C% 2 (y) (4.19) 



and 



quK-M = k~ M 1/2 f^d y cT{y)cT{y)c¥{y)c%\y) (4.20) 

arise from the expansion of the 2 and 3 terms, and the normalization factor ku is defined 
by Eq. (4.16). Note that Pu-m vanishes if the sum of any two indices is greater than the 
third, and that qNNN-,N > for all N. 
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Each term of the sum in Eq. (4.18) must vanish separately. For |A| sufficiently small, 
we expect to find a small amplitude solution in which a single coefficient, A N , is dominant. 
Explicitly, 



= AA N 
and 



PNN;NA N 



2b 



Pmn-,nA m A n - sgn (X)q NNN - N A 



|A| M^N 



(4.21) 



0=\J3- M(M + 3)] A M - -f=p NN . M A 2 N + 

A 



N, 



(4.22) 



where in each equation the omitted terms are higher order in A^. Using the second equation 
to eliminate A M from the first leads to 



= A 



Pnn; N A n - cA 2 N + 0(A%) , 



(4.23) 



where 



2b 2 2N 

C = Sgn (\)q N NN;N + TT7 H 



P-M(M + Z) 



PNN;MPNM;N 



(4.24) 



This gives 



1 

2~c 



bp 



NN;N 



± 



b 2 (p 



NN:N 



A 



+ 4cA 



(4.25) 



If the cubic term in Eq. (4.8) vanishes, 6 = 0, there are real solutions for A^ only if 
cA > 0, with An = ±^J A/c. Thus, for A > there are no small amplitude bounces if 
P < N(N + 3). As (3 is increased past this critical value, two solutions appear. For odd iV 
these are a physically equivalent true-to-false and false-to-true pair, while for even iV they 
are two distinct solutions, one true-to-true and the other false-to-false. Things are similar 
if A < 0, except that the solutions only exist for (3 < N(N + 3). 

For nonzero b and odd N, the situation is almost the same. Because Pnn-,n vanishes for 
odd N, the condition for existence of a solution is again that cA > 0, although now c and 
A may have different signs. As before, there are two distinct but physically equivalent small 
amplitude solutions, and these appear at the critical value j3 = N(N + 3). 

For nonzero b and even N, the situation is more complicated. From Eq. (4.25), we see 
that the critical value of (3 where solutions first appear is not (3 = N(N + 3), but rather 



13 = N(N + 3) 



4c|A 



(PNN;nY 



(4.26) 



which is a bit smaller (larger) if c is positive (negative). In contrast with the previous 
cases, when the solutions first appear they have a finite amplitude and both have the same 
sign for A N . As (3 is increased (decreased) through A^(A^ + 3), one of the solutions passes 
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through zero and then takes on the opposite sign. After this point the two solutions remain 
of opposite sign, with the magnitude of increasing for both. 

Some words of caution are in order here. Our analysis is valid only for sufficiently small 
An. If b/\j\\\ is too large, this condition may not be met when (3 is given by Eq. (4.26), 
in which case it will never be satisfied by the larger of the two solutions. However, as 
A — > 0, the other solution for A N tends to zero, and so should be reliable. As discussed at 
the beginning of this section, the zeros of h(<po) appear or disappear in pairs. Hence, the 
existence of one small amplitude solution implies the existence of the second. We therefore 
expect the qualitative features — the emergence of two solutions, with finite amplitude, at 
a value of (3 not corresponding to an integer — to be valid even if the quantitative analysis 
is not completely reliable. 

The results of this section can be compared with the analysis of Jensen and Steinhardt 
[14]. They argued that if \V"\ is monotonically decreasing as one moves from the top of the 
potential toward the false vacuum there will be a unique k = 1 Coleman-De Luccia bounce 
if (in our language) f3 > 4, and that this solution merges with the Hawking-Moss solution as 
P — > 4. When applied to Eq. (4.8), their condition on V" translates into the statement that 
A > 0. If b 2 /\ is not too large, c will also be positive, and our results will agree with theirs. 
However, for b 2 /X sufficiently large, c can be negative even though A is not. In this case, 
(3 = 4 is an upper bound for the existence of these solutions, and our results will disagree 
with their analysis. We believe that this can be attributed to the fact that their claims are 
based on the argument that the potential gets "flatter" as one moves away from the top of 
the barrier; for large b this is true, in the relevant range, only on one side of the potential. 

Perturbative methods similar to those used in this section have recently been applied to 
the calculation of the bounce action [15]. 

V. NUMERICAL RESULTS FOR A > 

Once the amplitude become large, the analysis based on a single dominant term in 
Eq. (4.17) fails, and we must resort to numerical methods. Before describing our numerical 
results, we outline some qualitative arguments that may provide some useful insight. We 
focus here on the case A > 0, where the quartic term in the potential is positive. We will 
return to the A < case in Sec. VI. 

We have seen that h(4> ) has a zero at (f) top , corresponding to the Hawking-Moss solution, 
for all values of (3. As (3 increases, pairs of new zeros, corresponding to bounce solutions 
with ever increasing values of k, appear at <f) top as (3 passes through successive critical values. 
When the analysis of Sec. IV C can be applied, these zeros move apart as (3 increases. We 
expect this behavior to continue even after the zeros have moved beyond the small amplitude 
regime. Thus, at any given (3 there should be a sequence of bounce solutions. For 6 = and 
N(N + 3) < (3 < (N + 1)(N + 4), there would be two solutions each for k = 1, 2, . . . , N, 
with the two solutions for each odd k simply being "y-reversed" versions of each other. For 
nonzero b the picture would be similar, except that the new solutions for even N appear 
somewhat before (3 = N(N + 3). Figure 3 gives a schematic illustration of the expected 
pattern of zeros. 

We can also make some predictions about the form of these solutions. It seems reasonable 
to expect that, even after <f> has moved well beyond the small amplitude regime, the average 
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FIG. 3. A schematic plot of the starting points of bounce solutions to a A > theory, such as 
that discussed in Sec. V, with (3 slightly greater than 28. The numbers next to the points represent 
the number of times the solution crosses the top of the potential barrier. 

period of the oscillations will continue to be determined to leading order by (3. Increasing 
(3 will decrease this period, and so should cause a plot of the oscillations of <fi to shrink 
horizontally at the same time that it is expanding vertically. With the number of oscillations 
in a given solution remaining fixed, this will open up "gaps" at the ends of the interval. In 
order to satisfy the boundary conditions that d(f)/dy = at y = and y = n, the field must 
be relatively constant in these gaps. This, in turn, implies that the field near the end points 
must be close to one or the other of the minima of the potential. As (3 is increased further, 
the fraction of the ^/-interval occupied by oscillations will continue to decrease, while the 
regions occupied by these false and true vacuum regions will expand. 

We have tested these ideas numerically by considering a theory governed by a quartic 
potential than can be written in the form 



where if) — <f>/v is a dimensionless field that has been rescaled so that for g = the minima 
of the potential are at ip — ±1; in the notation of Eq. (4.8), this corresponds to b = —g/3/2v 
and A = f3/v 2 . We focus here on our results for a single value of (3. However, we have 
explored solutions for a much wider range of values, going as high as 10 4 . In all cases, our 
results are consistent with the picture we describe here. 

In Fig. 4 we show the bounce solutions for g = \j2\[2 and (3 = 70.03, corresponding to 
N = 7. As expected, we find solutions for every integral value of k up to k = 7, with two 
distinct solutions for each even value of k. Although it is barely apparent on the figures 
(except for k — 7), it should be noted that, among solutions starting on a given side of the 
barrier, the initial value O monotonically approaches <p top as k increases, in accord with our 
expectations. 

The k — 1 solution is the Coleman-De Luccia bounce. For k = 2, there are two solutions. 
One starts and ends near the true vacuum, and might be viewed as a two-bounce solution 
analogous to the flat space multibounce solutions that are encountered in the path integral 
treatment. Indeed, its field profile is very close to that expected from two independent 
bounces and its action is approximately twice that of the single bounce. The other k = 2 
solution starts and ends near the false vacuum. It, too, might be interpreted as a two- 
bounce solution, if one took the Coleman-De Luccia bounce to represent decay from the 
true vacuum to the false. From this viewpoint, the false vacuum region corresponds to the 
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FIG. 4. Bounce solutions for the asymmetric (g = potential discussed in the text. For 

these solutions (3 = 70.03. 
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bubble interior; since this is large compared to the size of the four-sphere, there must be 
significant distortion to be able to fit in two bounces. Note that even in the region separating 
the two "interiors" does not get very close to its true vacuum value. 

The solutions with k > 2 all have intermediate oscillations about top - As expected, the 
frequency of these oscillations seems to be determined primarily by (3, with the separation 
between zeros of the field being roughly independent of k. The existence of intermediate 
oscillations seems to have little effect on the true vacuum bubble region, with the field 
profiles near y = n for the k < 7 false-to-true solutions and near y = and y = n for 
all the true-to-true solutions being very similar to that for the corresponding region of the 
Coleman-De Luccia solution. This is not so for the false vacuum regions. 

The plot for the k = 7 solution shows that the linear approximation works quite well, 
as should be expected with (3 so close to the critical value of 70. The numerical solution 
is indistinguishable to the eye from the Gegenbauer polynomial, and the initial value, <p , 
agrees well with the prediction of Eq. (4.25). 

In Fig. 5, we show the analogous results for the case of g — 0, where the two minima are 
degenerate. In flat spacetime, there would be no bounces for this case but, as we see, this 
is no longer so when gravitational effects are included. 

Finally, in Figs. 6 and 7 we show the actions for these solutions. (More precisely, what we 
actually plot in these figures is B = Sbouncc — SW, which is the quantity that would appear in 
the exponent if these solutions were interpreted as contributions to the process of tunneling 
out of the false vacuum.) As k increases, B approaches the Hawking- Moss value £>hm- This 
can be understood by noting that in the linearized approximation the contribution from the 
gradient terms exactly cancels that from the oscillatory part of the potential term, so that 
the total action is exactly equal to that of the Hawking-Moss solution. 

VI. BOUNCES WITH FLAT POTENTIAL BARRIERS 

Our analysis of small amplitude bounces in the previous two sections showed that for 
a wide class of potentials — those with a positive fourth derivative, and not too large a 
third derivative, at top — there is neither a small oscillation bounce nor a bounce smoothly 
related to a small oscillation bounce if | V /r (^>top) | < 4H 2 . This result is intuitively quite 
plausible. One might expect | V" (4> tov )\~ l l 2 to set a natural scale for the size of the bounce, 
so that if this were too large, relative to the radius H' 1 , the bounce would not fit on the 
four-sphere. 

However, our results for the case where the fourth derivative (or more precisely, the 
related quantity c) is negative show that this seemingly plausible argument cannot be quite 
correct. For this case we saw that there was a small amplitude k = 1 bounce only for 
a range of values extending downward from (3 = 4. One logical possibility is that these 
solutions continue (although not necessarily with small amplitudes) all the way down to 
(3 = 0. Alternatively, there might be some minimum value (3 m i n (which could depend on 
the other parameters of the theory) for which the bounce exists. Since solutions can only 
appear and disappear in pairs, this latter possibility would require the existence of a second 
solution, which would also appear at /5 m i n but which would persist beyond (3 = 4. 

In this section we will explore this direction in more detail. Because the potentials that 
evade the (3 > 4 bound are characterized by being relatively flat at the top of the barrier, 
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FIG. 5. Bounce solutions for the symmetric (g = 0) potential with two degenerate vacua that 
is discussed in the text. Again, (3 = 70.03. 
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FIG. 6. The tunneling exponent B = -S^bounce) — S(4>{ v ) for the asymmetric potential bounce 
solutions shown in Fig. 4. The values are given as fractions of the Hawking-Moss value, and are 
plotted as a function of the number of times the solution crosses the top of the barrier. Boxes 
correspond to solutions which interpolate between the true and false vacuum sides of the barrier, 
diamonds correspond to "true-to-true" solutions, and triangles to "false-to-false" ones. 
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FIG. 7. The tunneling exponent B for the symmetric potential bounce solutions shown in 
Fig. 5. The values are given as fractions of the Hawking-Moss value, and are plotted as a function 
of the number of times the solution crosses the top of the barrier. Boxes correspond to solutions 
which interpolate between the the two vacua, and triangles correspond to solutions which begin 
and end on the same side of the barrier. 
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we begin by studying a toy model, denned by the potential 

f CO + a) + V , <p < -a 

V = 



V , —a < < a 
{-C(<f)-a)+V , 



(6.1) 



> a 



that is absolutely flat at the top. This potential has the additional advantage that the field 
equation can be solved analytically. 4 

Inserting this potential into Eq. (4.4), and defining 



leads to 



dx 
dy 



= < 



X{y) = sin y— , 

' 07 sin 3 y , < —a 
, —a < < a 



(6.2) 



(6.3) 



-07 sin y . 



> a 



where 



7 



C 



(6.4) 



Because C/a gives a rough measure of the average value of \V"\ for our toy potential, 7 is 
somewhat analogous to the quantity (5 that appeared in the analysis of the two previous 
sections. 

We first seek a k — 1 (Coleman-De Luccia) bounce. The symmetry of the potential 
implies that this solution (as well as any other bounce with odd k) must have 0(7r/2) = 0. 
Let 0(0) = 0o < —a, and define y\ < n/2 by the condition <f>(yi) = —a. Integrating Eq. (6.3) 
over the interval < y < y\ gives 



x(yi) — a l I dusm u. 
Jo 



The fact that x is constant over the flat part of the potential implies that 

dv 



a = 0(tt/2) - 0( yi ) = X (yi) f 



i yi sin v 

Combining these two equations gives the requirement that 



(6.5) 



(6.6) 



4 The fact that this potential has no minima is irrelevant for our purposes, since the bounce 
solutions never quite reach either the true or false vacuum in any case. One could, of course, 
modify the potential to give minima at large values of |0|. 
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0.25 0.5 0.75 1 1.25 1.5 

FIG. 8. The functions iq(yi) and F 2 (yi) for the toy model defined by Eq. (6.1) 




FIG. 9. The two k = 1 bounces for our linear toy model when 7 = 10. 
1 rvi /" 7r /2 dy 

-=/ dusin 3 u — »- = Fi(«i). (6.7) 

Once a value of yi satisfying this condition has been found, it is a straightforward matter to 
find 0o by integrating back from y\ to y = and to then integrate forward to obtain (f>(y) 
for all y. 

From the plot of F\{y) shown in Fig. 8, we see that Eq. (6.7) has no solutions, and hence 
there is no Coleman-De Luccia bounce, if 7 < 1/Max [-Fi(y)] ~ 7.70. For any 7 greater than 
this critical value there are two acceptable values for y 1; and thus two distinct Coleman-De 
Luccia bounces. As examples of this, in Fig. 9 we show the two bounces solutions for 7 = 10. 

In general, these values for y\ can only be found numerically. However, approximate 
expressions can be obtained for the limiting case 7 ^> 1. In one solution, with 
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Vi ~ \ - 0o~-a(l + -], (6.8) 

the field is mostly on the flat part of the potential. In the other solution, with 

3 , /l + ln4\ nnn ,„ n . 

yi « - - — 0o ~ - I g I 07 « -.398a7 , (6.9) 

the field is mostly on the sloping parts of the potential. The latter solution has the lower 
action. 

For k = 2 (or any even value of k), d(f>/dy vanishes at y = tt/2. Let us define 2/1 < y 2 < 
tt/2 by 0(yi) = —0(?/2) = — a- Integrating Eq. (6.3) then gives 

x(j/i) = 07 / cousin -u (6.10) 
•/ 



/■tt/2 

xivz) = a l dusin 3 u. (6.11) 



Because the field is in the flat part of the potential for y\ < y < y 2 , these integrals must be 
equal; this implicitly defines y 2 as a function of 2/1. 
By analogy with Eq. (6.6), we also have 



ry2{yi) du 

2a = 0(y 2 ) - 0(^/0 = / = 2a~fF 2 ( yi ) . (6.12) 

»/ yi Sill W 



Referring again to Fig. 8, we see that there are no solutions for 7 < 1/Max [-^2(2/)] ~ 21.7, 
but two solutions for all larger values of 7. In the cases where there are solutions, these have 
values of 2/1 (and also of O ) that fall between the corresponding k = 1 values. In the limit 
7>1, the k = 2 solutions have 

j/! ss J— 0o~-a(l + -| (6.13) 

V 7 V 7/ 

and 



2/1 w 1.22 0o w -.319o7 . (6.14) 

This procedure can be extended to all higher values of k. For any given k, solutions only 
exist if 7 is larger than a critical value 7 cr (/c); if 7 > 7 cr (^) there are two allowed values for 
0o, both of which fall between the allowed valued for the k — 1 solutions with the same 7. 
For k — 3, the critical value is 42.7. More generally, 7 cr (^) is an increasing function of k and 
is 0(k 2 ) for large /c; this behavior is reminiscent of the critical values of (3 that govern the 
appearance of small amplitude solutions in potentials with a positive quartic term. 

The situation is illustrated in Fig. 10, which gives a schematic summary of the possible 
values for O for a given value of 7. As 7 is increased, the zeros marked by closed circles move 
outward, while those marked by open circles move inward. When 7 reaches a critical value, 
two new pairs of coincident zeros (one pair on each side of the potential) appear and then 
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FIG. 10. A schematic graph of the starting points of bounce solutions to our linear toy model. 
Closed circles represent solutions which move outward as H is decreased, and open circles represent 
solutions which move inward. The numbers next to the points represent the number of times the 
solution crosses the top of the potential barrier. 



begin to separate. With increasing 7, accumulations of open circles develop near <p = ±a, 
while the outwardly moving closed circles become further and further apart. 

The actions of these solutions may be either greater or less than that of the Hawking- 
Moss solution. For example, the k — 1 solution with the smaller |0o| always has a higher 
action than both the other k = 1 solution and the Hawking-Moss. The action of the other 
k — 1 solution is greater than Shm for small 7, but not for large 7. 

This toy model has the advantages of allowing analytic solution and of emphasizing 
that it is possible to have a Coleman-De Luccia bounce with an arbitrarily small value of 
|V"(0top)|- However, one consequence of its somewhat artificial form is that the bounce 
solutions have some special properties that we would not expect in a more generic potential. 
For example, the values of O that move inward with increasing 7 (i.e., those indicated by 
open circles in Fig. 10) accumulate near the "corners" of the potential at ±a, but never 
move further inward where they might meet and "annihilate" . 

A less artificial example can be obtained by considering a potential of the form V = 
—a<f) 2 — bcj) A + c0 6 , with a, b, and c all positive. With three parameters to vary, it is possible 
to make the potential relatively flat, as defined by a scale-independent criterion, at the top 
of the barrier. For example, if ±v are the minima of the potential, we can define an averaged 
second derivative V^ g by 

Cs = yW2)-y ( -,/2) (6i5) 

By choosing b 2 ^> ac we can make |V"(0)| <C |V^ g |. 

We have followed the evolution of the bounce solutions for such a potential (with b 2 /ac = 
80) as H is varied with the parameters in the potential held fixed. For very large H, there are 
no bounce solutions, in analogy with the low (3 or low 7 behavior of our previous examples. 
As H is decreased, two pairs of zeros, corresponding to initial values for k — 1 bounces, 
appear at nonzero points on opposite sides of the barrier, just as in our toy model. (In the 
example we studied, these first appeared at \V^ g \/H 2 = 4.03.) As H is reduced further, 
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FIG. 11. A schematic graph of the starting points of bounce solutions for the sixth-order 
potential discussed in the text. Closed circles represent solutions which move towards one of the 
vacua as H is decreased, and open circles represent solutions which move towards the top. The 
numbers next to the points represent the number of times the solution crosses the top of the 
potential barrier. 

each pair splits, with one zero moving inward and one outward. 5 With further decreases in 
H, additional zeros, corresponding to bounces with k = 2, k = 3, etc. appear, just as in the 
toy model. 

In contrast with the toy model, the inwardly moving zeros proceed all the way in to 
= 0, reaching that point exactly when predicted by the analysis of Sec. IV C. Thus, 
two k = 1 zeros meet and annihilate when H 2 = a/2, two k = 2 zeros annihilate when 
H 2 = a/5, etc. We have also encountered a more complex type of behavior for some values 
of the parameters. As described above, two pairs of k = 2 zeros appear at a certain critical 
value of H. After moving inward some distance, each of the inner zeros splits into three 
zeros. These then separate for a time, and then rejoin to form a single zero, which eventually 
moves to = 0. A schematic snapshot of this example for a fixed value of H is shown in 
Fig. 11. 

The examples considered in this section clearly demonstrate that the existence of a 
Coleman-De Luccia bounce need not impose a lower bound on |V"(0top)|- However, they sug- 
gest that a more generalized version of this bound, involving an averaged value of \V"\/H 2 , 



5 Another example whose k = 1 solutions display a similar behavior was discussed by Jensen and 
Steinhardt [16]. 
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does hold. This is consistent with the results of [17], where it is argued that a necessary 
condition for the existence of a Coleman-De Luccia bounce is that \V"((p)\/H 2 > 4 for some 
value of 0. 

VII. CONCLUDING REMARKS 

In this concluding section we will first summarize our results and then discuss how they 
might be affected if we relax some of the restrictions that we have imposed to simplify the 
analysis. We will then briefly discuss the physical implications and interpretation of these 
solutions. 

A. Summary 

In this paper we have studied in some detail a class of oscillating bounce solutions that 
arise when the effects of gravity on vacuum tunneling processes are taken into account. 
These are qualitatively distinct from the Coleman-De Luccia and Hawking-Moss solutions. 
In fact, they can be viewed as, in a sense, interpolating between these two limiting cases. 

These oscillating bounces exist on a space that is topologically a Euclidean four-sphere. 
They are 0(4)-symmetric, so that the scalar field depends only on the parameter £, which is 
proportional to the azimuthal angle. At k > 1 values of £ the scalar field takes on its value 
at the top of the potential barrier. These values correspond to k three-spheres that can be 
viewed as dividing the four-sphere into two end caps and k — 1 bands. Typically, in each of 
the end cap regions the field approaches close to, but does not quite reach, one of the two 
possible vacuum values. In the intermediate band regions the field is alternately on the true 
or false vacuum side of the barrier, although it does not generally wander very far from the 
top of the barrier, with this deviation decreasing as k increases. 

For k — 1, the solution is simply the Coleman-De Luccia bounce. The k = 2 solution for 
which the field in both end caps is in the true vacuum region can be viewed as a two-bounce 
solution. However, the solutions with higher values of k are not in any sense multi-bounce 
solutions. For any fixed value of the parameters of the theory, there are only a finite number 
(possibly zero) of these solutions. When solutions exist, they have values of k taking on 
all integer values from 1 to some k max , with possibly several distinct solutions for the same 
value of k. The determination of precisely how many solutions exist is closely related to the 
older problem of determining when a potential admits a Coleman-De Luccia bounce. 

In a rough sense, the controlling quantity here is the ratio of the second derivative of 
the potential to the value of H 2 . This is physically plausible, since for "typical" potentials 
V" is related to the characteristic mass scale, which in turn determines the characteristic 
size of a classical solution. If this last quantity were too large, one might imagine that the 
bounce would not "fit" on the four-sphere. Early discussions focused on (3 = \V"((j) t0 p)\/H 2 : 
and argued that there was no Coleman-De Luccia bounce if f3 < 4. We have shown that 
for a broad class potentials (5 is indeed the controlling parameter and that /c max = N, where 
N is the largest integer such that N(N + 3) < /3. As j3 is increased through one of these 
critical values, two (possibly equivalent) k = N solutions appear. Precisely at the critical 
value, the new solutions are identical to the Hawking-Moss solution, and so can be thought 
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of as being "all wall" . (As described in Sec. IV, matters can be slightly more complicated if 
k is an even number.) 

Not all potentials give this behavior, as is evident from our results for the case where 
the fourth derivative of the potential is negative at the top of the barrier. Here, the critical 
values of f3 give upper, rather than lower, limits for the existence of small amplitude bounce 
solutions. The situation is somewhat clarified by considering potentials where the barrier 
is atypically flat. For these, the parameter whose value determines the first appearance of 
a bounce is not the second derivative at the top of the barrier. Instead, one can use an 
averaged value of \V"\/H 2 (whose precise definition may vary from case to case) to define a 
quantity 7. As 7 is increased through a /c-dependent critical value 7 1; a pair of Absolutions 
appears. These coincide, and are nontrivial, when they first appear. In the simplest case, 
one of these will persist for all 7 > 71, while the other survives only for a finite range 
71 < 7 < 72, with the parameters that give 7 = 72 also giving (5 = k(k + 3). When this 
critical point is reached, the latter solution and its ^-reversed image annihilate and merge 
into the Hawking-Moss solution. More complicated patterns are also possible, but must 
satisfy the constraint that solutions only appear and disappear in pairs. To satisfy this 
requirement there must be at least one solution that survives for arbitrarily large 7, and a 
second solution that either merges with the Hawking-Moss at a finite value of 7 or, as with 
our toy model, also persists for arbitrarily large 7. 

B. Robustness of our results 

To simplify our analysis, we have imposed two restrictions. First, we have only considered 
solutions with 0(4) symmetry. Second, we have assumed that the potential is such that we 
can use the fixed background approximation. We must now ask how robust our results are, 
and how they might change were these restrictions relaxed. 

The restriction to 0(4)-symmetric configurations is the common practice in flat space 
treatments of vacuum tunneling. For single-field potentials satisfying a rather plausible set of 
conditions, it can even be shown that the lowest action bounce has 0(4) symmetry. Indeed, 
although there are nonsymmetric approximate stationary points corresponding to several 
widely separated bounces, it may well be that most theories admit no bounce solutions 
without 0(4) symmetry. 

This situation is no longer the case when gravitational corrections are included, largely 
because the Euclidean space is now a compact manifold. To see this, consider the regime 
where the true vacuum bubble of the Coleman-De Luccia bounce has a radius much less than 
if -1 . Now arrange several such bubbles in a maximally symmetric fashion around the four- 
sphere (e.g., at the six vertices of a five-dimensional hypertetrahedron) . It is clear that some 
mild local smoothing of the configuration will give a stationary point of the action, and hence 
a bounce solution without 0(4) symmetry. The solution will, however, be invariant under a 
discrete subgroup of 0(5). It is quite possible that there might be additional solutions, with 
oscillatory behavior, with the same discrete symmetry. However, there no reason to expect 
more than a finite number of these. Further, we would expect that the actual number of 
such solutions would have a dependence on the parameters of the theory that is similar to 
that found for the 0(4)-symmetric solutions. 
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A more significant limitation on our analysis is the use of the fixed background approxi- 
mation, which is clearly inapplicable to many potentials of interest. Nevertheless, we expect 
most if not all of the qualitative, and indeed many of the quantitative, features of our re- 
sults to persist. Thus, for solutions that begin as small deviations from the Hawking- Moss 
solution when (3 reaches a critical value (either from above or below), the value of (3 cr will 
be exactly as in our analysis, provided that H is calculated using V(<f>to P )- 

A second case to consider is the toy model of Sec. VI. Let us assume that the linear falloff 
of the potential is terminated at some finite value of \<f>\ in such a way that V((f>) remains 
everywhere positive. Further, let 7 still be defined using the value of the potential at the 
top of the barrier. The fact that V{4>) is actually lower than this in some regions would have 
the effect of decreasing the effective value of H; from Eq. (3.4), we see that the nonzero <fi 
acts in the same direction. The Euclidean space will therefore be larger, making it easier 
to fit in a bounce solution. As a result, we expect a reduction in the values of 7 at which 
solutions for a given k will first appear, and thus an even larger range of parameters that 
allow a bounce with an absolutely flat potential. 



C. Physical interpretation 

We now turn to the physical interpretation of these solutions. As with the Coleman-De 
Luccia and Hawking-Moss bounces, the hypersurface passing through £ = and £ = £ max 
gives the initial data for the classical evolution after the vacuum transition governed by 
the bounce. This hypersurface is a three-sphere of radius ~ H -1 , apparently implying 
that the transition takes place at the "waist" of the de Sitter hyperboloid. This is generally 
understood to be an artifact of the formalism, with the physical results instead being adapted 
to a single horizon volume. 

Formally, the classical evolution can be obtained by an analytic continuation of the 
bounce solution. In doing so, particular care must be taken at the points £ = and £ = £ max 
of the Euclidean solution. After rotation to a Lorentzian signature these become the lightlike 
boundaries of two antipodal regions on the de Sitter hyperboloid [13]; this generalizes the 
manner in which the origin of the flat space Euclidean solution becomes a lightcone after 
the Wick rotation to Minkowski space. 

In the case of the Coleman-De Luccia bounce, this procedure yields a Lorentzian space- 
time in which portions of two de Sitter spaces, corresponding to the two vacua, are separated 
by a well-defined bubble wall. When the initial bubble size is small, the bounce can be viewed 
as having an initial state hypersurface corresponding to the false vacuum (see Fig. 2), and the 
vacuum transition is most naturally interpreted as quantum mechanical tunneling. When 
the initial bubble size is comparable to H -1 , the initial state is not evident in the bounce, 
but is instead determined only by which vacuum action is subtracted in the calculation of 
B. The transition thus has a more thermal character, reflecting the existence of a nonzero 
de Sitter temperature T dS . 

The Hawking-Moss solution is even more clearly thermal in character. Here the initial 
data for the classical evolution is a homogeneous region with the field balanced at the top of 
the potential barrier. Although, strictly speaking, the classical Lorentzian evolution would 
leave the scalar field at the top of the barrier forever, this solution is unstable against small 
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fluctuations and so will break up, in a stochastic fashion, into regions that evolve toward 
one or the other of the two vacua. 

The oscillating bounces yield a hybrid of these two cases. The end cap regions near £ = 
and £ = £ max evolve into two vacuum regions bounded by well-defined walls; these vacua 
can be the same or different, and have no definite relation with the state before the vacuum 
transition. In the intermediate region, the scalar field profile oscillates about the top of the 
barrier 6 . With initial conditions given precisely by the bounce solution, these oscillations 
would be preserved under the classical evolution, giving a large and exponentially expanding 
region in which the field is near the top of the barrier. However, just as with the Hawking- 
Moss solution, there will be instabilities against small fluctuations. We expect these to 
lead to a breakup into a stochastic mixture of regions evolving toward the true and false 
vacua, separated by relatively narrow transition regions. (We differ here from the scenario 
advocated in Ref. [6].) 

The relative importance of these three types of solutions depends on the parameters of 
the theory. For "typical" potentials, such as that studied in Sec. V, the various regimes can 
be characterized by the value of j3. (If the potential is atypical, e.g., by being unusually flat at 
the top of the barrier or by having two very different mass scales, the regimes will be defined 
by a more complex combination of parameters, but the discussion would be analogous.) 
When j3 >> 1, there will be a Coleman-De Luccia bounce, a Hawking- Moss solution, and 
many oscillating bounces. However, the actions of the Hawking-Moss solutions and the 
k > 2 oscillating bounces will be many times greater that of the Coleman-De Luccia bounce, 
which will clearly dominate. This is a regime of quantum tunneling transitions followed by 
deterministic classical evolution. At the other extreme is the case /3 ^ 1, where only the 
Hawking-Moss solution is the only bounce. This is a regime of thermal transitions followed 
by stochastic evolution. Intermediate between these is a regime of moderate (5 that admits 
the Hawking-Moss, one (or more) Coleman-De Luccia, and several oscillating bounces. It is 
in this transitional regime that the oscillating bounces are most likely to play a role. 

Finally, we should comment on the question of negative eigenmodes. In the case of flat 
space quantum tunneling, the existence of a single negative mode (or, possibly, a number 
equal to one mod four) was an essential ingredient to the path integral derivation of Y. 
Indeed, it can be shown [18] that for a wide class of theories in flat space the bounce of lowest 
action has only one negative mode. The situation in curved space, and in particular for the 
oscillating bounces, needs further clarification, both as to how these affect the contribution 
to the vacuum transition rate and to their role in the stochastic breakup of the intermediate 
regime after the transition. We plan to return to this subject in a later publication. 
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6 An exception is the k = 2 "true-to-true" solution, in which the intermediate region is approxi- 
mately false vacuum at nucleation. 
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